IFUP-TH 29/93 



Two dimensional SU {N) x SU (N) chiral models on the lattice 



Paolo Rossi and Ettore Vicari 
Dipartimento di Fisica deU'Universita and I.N.F.N., 1-56126 Pisa, Italy 



Q"^ , Lattice SU{N) x SU (N) chiral models are analyzed by strong and weak coupling 

expansions and by numerical simulations. 

12^^ order strong coupling series for the free and internal energy are obtained 

for all > 6. Three loop contributions to the internal energy and to the lattice 

■ /^-function are evaluated for all N and non-universal corrections to the asymptotic A 

OO . parameter are computed in the "temperature" and the "energy" scheme. Numerical 

I simulations confirm a faster approach to asymptopia of the energy scheme. A phe- 

^ ' nomenological correlation between the peak in the specific heat and the dip of the 

' /3-function is observed. 

{•"•^ ■ Tests of scaling are performed for various physical quantities, finding substantial 

, scaling at ^ ^ 2. In particular, at = 6 three different mass ratios are determined 

[ numerically and found in agreement, within statistical errors of about 1%, with the 

0^ I theoretical predictions from the exact S-matrix theory. 
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I. INTRODUCTION 



Two dimensional SU{N) x SU{N) principal cliiral models defined by the continuum 
lagrangian 

L = ^Trd.Ud.U^ (1) 

are the simplest quantum field theories sharing with non-abelian gauge theories the property 
of asymptotic freedom and whose large N limit is a sum over planar diagrams. Due to the 
existence of higher-order conservation laws, multiparticle amplitudes are factorized, and 
exact S matrices have been proposed The resulting bound state mass spectrum is 

represented by 

M^^M'j^^ . 1<.<^-1 , (2) 

and the bound state of r particles transforms as the totally antisymmetric tensor of rank r. 
The mass-A parameter ratio has been computed, and the result is Q 

M /Stt sin tt/N , , 

- ' • (3) 



his V e vr/iV 

A "standard" lattice version of principal chiral models is obtained by choosing the action 
S = -2/3A^^ReTr[[/.f/I+^] , /? = -L , (4) 

n,iJ. 

whose properties have been investigated by several authors ||5 Hll|] especially by strong- 
coupling and mean field methods. Numerical simulations have also been performed (at 
= 3), most recently by Dagotto-Kogut [1^ and Hasenbusch- Meyer |[T3[ . 



As a preliminary step within a more general program whose ultimate goal is performing 
the numerical 1/A^ expansion of matrix- valued field theories, we decided to explore the 
properties of principal chiral models at larger-than-usual values of A^. In particular we 
wanted to investigate the following issues: 

- region of applicability, accuracy and A^ dependence of the strong coupling series. 

- onset of scaling, with special attention to the interplay between thermodynamical (peak 
in the specific heat) and field theoretical (dip in the /3-function) effects. 

- check of conjectured exact results (especially mass ratios) by Monte Carlo measurements 
in the scaling region. 

- role of coupling redefinitions in the widening of the asymptotic scaling regions. 

To this purpose, we performed a variety of strong coupling and weak coupling calcula- 
tions, and a number of numerical simulations for different values of A^, and especially at 
A^ = 6, where the mass spectrum is sufficiently non-trivial (two independent mass ratios 
can be measured and compared with prediction), and 0(1/A^^) effects should be already 
significantly depressed. 

In the present paper we only report on our analytical results, without offering any details 
on the derivations, that will be presented elsewhere. 
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II. ANALYTICAL RESULTS 



A. Strong coupling 



We found that the most convenient approach to strong couphng is the character expan- 
sion. Free energy character expansion for U{N) chiral models to twelfth order and mass 
gap expansion to fifth order were presented in Ref. 0]. The formal extension of these se- 
ries to SU{N) is easily achieved with the abovementioned precision for N > 6. Paying 
some attention in order to avoid double-counting, SU (6) can be also obtained by the same 
technique. 

We found explicit representations of the coefficients of the SU{N) character expansion 
in the strong coupling regime in terms of Bessel functions, by generalizing the technique 



discussed in Refs. 



]. These representations are exact up to O 



2N 



. As consequence 



we could compute the SU (N) free energy to twelfth order in P for > 6 in two dimensions 



1 
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In the case = 6 an analysis of the O \j3^^^j and O 
1 „ 



'2N 



contributions led to the result 



2(5'' + 2(3" + 11.2/3'' + 68.602449^* + 374.945306/3^° + 6395.760105/3^^ + 



?12 



iV2- 



(6) 



The internal energy (per link) density E is immediately obtained from the previous results 
by 



E 



1 dF 



(7) 



These results have been used to draw the strong coupling curves in our figures and compare 
very well with numerical simulations in the region (3 ^ 0.25. 



B. Weak coupling 

Short weak coupling series for the free-energy density of U {N) and SU{N) chiral models 
were presented in Ref. 0. 

We calculated the energy density up to three loops finding 
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1 N'^ — 1 
E = I - ( —Re Tr \UnUL^ ) = — 



(3 (3^ 



(8) 



where 



_N^-2 



3iV^ - 14iV2 _^ 20 N^-4:N^ + 12 N^-8N^ + 24: 

^^^^^ ^ 64iV4 ^ 



'2 , 



(9) 



"2 768^r4 + 64JV4 + 64iV4 ^ 

Qi and Q2 being numerical constants: Qi = 0.0958876 and Q2 = —0.0670. 

Asymptotic scaling requires the ratio of any dimensional quantity to the appropriate 
power of the two loop lattice scale 



K21 = (6oT)-^^/^«exp(-^ 



(10) 

to go to a constant as T — 0. 60 and bi are the first universal coefficients of the expansion 
of the /3-function: 

= -a^T = -boT'' -biT^ - 62, + ... (11) 



- 8^ ' - 128^^ 



(12) 



Evaluation of the ratios of A parameters requires a one loop calculation in perturbation 
theory, which leads to M: 



At. 



32 exp TT- 



N^-2 ' 
2m 



(13) 



In order to get a more accurate description of the approach to asymptotic scaling we 
performed the change of variables suggested by Parisi |jl5|, defining a new temperature Te 
proportional to the energy: 



IN 



- 1 



E 



NTf 



(14) 



Notice that the corresponding specific heat is, by definition, constant. The ratio of A^;, the 
A parameter of the scheme, and A^ is easily obtained from the two loop term of the 
energy density: 



Ae ( N'^-2' 

- — = exp TT- 



4Ar2 



(15) 



We encountered the usual (and yet unexplained) phenomenon of a much better conver- 
gence to asymptotic scaling for quantities plotted as functions of Pe |T^-|r8|. We tried to 
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check for a perturbative explanation of this phenomenon by computing the first perturbative 
correction to the two loop lattice scale 



A 



1/2 ^-877/3 



ml 



(16) 



in the standard and the /Se scheme, which requires the calculation of the three-loop term of 
the /3-function in both schemes. 

In the standard scheme we found 



b2r 



(27r)3 128 



1 + vr- 



2N^ 



TT 



'2N^ - 13A^2 ^ 18 



+ 4Gi 



(17) 



where d = 0.04616363 p|,|Ol. The equivalence of the SU{2) x SU{2) chiral model to the 
0(4) a model allows a check of this equation, indeed for = 2 it must give (and indeed it 
does) the same of the standard lattice 0(4) a model |]19 . 



The /9-function of the Pe scheme can be written in the form 



PEiTE 



where 



da 



C{T) 



sm 



m - 1 



1 dE 



C{T)h{T) 



(18) 



(19) 



is the specific heat and T must be considered as a function of Te- Expanding perturbatively 
Eq. (|18D and using Eq. (||) one finds 



^2e 



+ N\ (as - af) + Nbiai 



(20) 



As one may easily verify, the linear corrections to the two loop lattice scale in Eq. (16) 
are small and of the same order of magnitude (although of opposite sign). They cannot 
therefore explain the failure of the first and the success of the second scheme with respect 
to achieving asymptotic scaling. We believe that the origin of this phenomenon is fully 
non-perturbative, and it can presumably be traced to the phenomenologically apparent 
correlation existing between the peak in the specific heat and the dip in the lattice f3- 
function: the non-perturbative variable transformation that flattens the peak manages to 
fill the dip, in a theoretically yet uncontrolled way. 



III. NUMERICAL RESULTS 

We performed Monte Carlo simulations of the lattice SU{N) x SU{N) chiral models for 
a wide range of values of A^ (in particular A^ = 3, 6, 9, 15) and (3. Summaries of the runs are 
presented in Tables |,|I|,|ITT| and |rV|. 



In our simulations we implemented the Cabibbo-Marinari algorithm |21| to upgrade 



SU (N) matrices by updating its SU (2) subgroups. In most cases, we chose to update the 
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A^— 1 diagonal subsequent SU{2) subgroups of each SU{N) matrix variable by employing the 
over-heat-bath algorithm (for the "heat bath" part of it we used the Kennedy-Pendleton 
algorithm [p^ ). 

An important class of observables of the SU{N) x SU{N) chiral models can be con- 
structed by considering the group invariant correlation function 

G{x -y) = { ^Re Tr [f/(x)[/(i/)t] ) . (21) 

We define the correlation function from the second moment of the correlation function 
G{x). On the lattice 



4sin^ tt/L 



G(0,1) 



(22) 



where G{kx, ky) is the Fourier transform of G{x). The inverse mass gap is extracted from 
the long distance behavior of the zero space momentum correlation function constructed with 
G{x). Moreover we measured the diagonal wall- wall correlation length to test rotation 
invariance. M = 1/C,w should reproduce in the continuum limit the mass of the fundamental 
state. The first definition of correlation length C,g offers the advantage of being directly 
measurable, while the calculation of requires a fit procedure. On the other hand, since 
is an off-shell quantity an analytical prediction exists only for the inverse mass-gap (Eq. (^). 

In Tables |,p|,pT| and |3 we present data for the energy density E, the specific heat 
C = jj^, the magnetic susceptibility Xm defined from the correlation function G{x), the cor- 
relation length ^G, the dimensionless ratios ^g/^w and ^d/iw, respectively for = 3, 6, 9, 15. 

We carefully checked for finite size effects. It turned out that for 2 = L/^^ ~ 8 the finite 
size systematic errors in evaluating infinite volume quantities should be safely smaller than 
1%, which is the typical statistical error of our data. 

In Figs. and ^ we show the energy density versus (3 respectively at A^ = 6 and A^ = 9. 
There the strong coupling series up to twelfth order in f3 and the weak coupling one up to 
third order in f]~^ are drawn. 

As in other asymptotically free models, at all values of A^ the specific heat shows a peak, 
connecting the two different asymptotic behaviors: monotonically increasing in the strong 
coupling region and decreasing at large /?. In Figs. and | C is plotted respectively for 
A^ = 3, 6, 9, 15 with the corresponding 13*^^ order strong coupling series (except for A^ = 3). 
Increasing A^, the peak moves slightly towards higher P values {Ppeak — 0.285 at A^ = 6, 
Ppeak — 0.30 at A^ = 15), and becomes more and more pronounced. We found the position of 
the peak to be more stable at large A^ when plotting C versus ^g, as in Fig. ^ Notice that, 
increasing A^, the specific heat around the peak does not show any apparent convergence to 
a finite value, which might be an indication of a (first order?) phase transition at A^ = 00. 

The 12*'' order (13*'' order) strong coupling series of the energy (specific heat) are in 
quantitative agreement (within our statistical errors) for j3 ^ 0.2, and in qualitative agree- 
ment up to the peak of the specific heat, whose position should give an estimate of the 
strong coupling convergence radius. 

Tests of scaling, based on the stability of dimensionless physical quantities (for example, 
the ratio ic/Cw) and rotation invariance (checking that Cw/^d — 1); showed that, within our 
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statistical errors, the scaling region is reached already at small correlation lengths, i.e. for 
~ 2. Fitting data in the scaling region to a constant we found 



0.987(2) 


for 


= 


3 


0.993(2) 


for 


= 


6 


0.995(3) 


for 


A^ = 


9 


0.994(4) 


for 


A^ = 


15 



(23) 

Notice that scaling is observed even before the peak of the specific heat. Since strong 
coupling series should be effective in this region, it might be possible to calculate continuum 
physical quantities by strong coupling techniques. In order to investigate this issue work to 
extend the strong coupling series is in progress. 

We checked asymptotic scaling according to the two loop formula (p!0D by analyzing 
Mc/Al^21 = 1/('CgAl,2z)- In Figs. and |^ we show the corresponding data respectively 
for A^ = 3, 6, 9, 15. At all values of A^ we observe the usual dip in the /3-function, which is, 
again, more and more pronounced when increasing A^. Since — we compare Mq/ Al,21 
directly with Eq. (||) (using also Eq. (|13D), whose predictions, 

M/Ai 



48.266... 


for 


A^ = 


3 


72.412... 


for 


A^ = 


6 


77.989... 


for 


A^ = 


9 


81.001... 


for 


A^ = 


15 



(24) 

are represented by dashed lines in the figures. Notice that the Monte Carlo data are much 
larger than the predicted values, while the first perturbative corrections in Eq. (|T6|) are, in 
all cases, about 20% at /3 ~ 0.3. Furthermore, data show a similarity with the behavior of 
the specific heat, strengthening the idea of a strong correlation between the two phenomena. 

The approach to asymptotic scaling gets an impressive improvement using the (3e scheme. 
In Figs. and ^ we also plot 

Mg/Kl,2i\e = — ^ X ^ . (25) 

Now data approach the correct value, and the discrepancies are even smaller than the linear 
correction calculated in Sec. fT| (which is about 15% at /3 ~ 0.3). So flattening the peak of 
the specific heat by performing the coupling redefinitions T ^ Tg, the dip of the /3-function 
disappears. We believe this to be the key point of the success of the Pe scheme in widening 
the asymptotic scaling region. The peak of the specific heat should be explicable in terms of 
complex /^-singularities of the partition function close to the real axis The sharpening 
of the peak with increasing A^ would indicate that the complex singularities get nearer and 
nearer to the real axis, pinching it at A^ = cxd where a phase transition is expected. Such 
singularities should also cause the abrupt departure from the weak coupling behavior. Then 
a coupling transformation eliminating the peak should move the complex /^-singularities 
away from the real axis, and therefore improve the approach to asymptotic scaling. 

From the Monte Carlo data and the exact result (^) we can extract the effec- 
tive A-parameters Al{N,P) and Ae{N, Pe)- Fig- §| and ^ show respectively the ratios 
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Kl{N,[3)/Kl,2i{N,I3) and Ke{N, Pe) / ^e,2i{N, [3e), where Kl,2i{N,P) and KeAN^Pe) are 
the corresponding two loop functions: Kl^2i{N^x) = Ae,2i{N,x) = (Svrx)"'^/^ exp(— Svrx). 
Similarly to the specific heat, the effective A-parameter AL{N,f3) does not give evidence of 
convergence at large N. On the contrary Ae{N,Pe) appears to approach a finite function 
Ae{oo, Pe), which is well approximated by the two loop formula. 

In conclusion, scaling and asymptotic scaling (in the Pe scheme) are observed at all 
values of considered, even around the peak of the specific heat. It is interesting to notice 
that, even though the behavior of the specific heat with respect to suggests the existence 
of a phase transition at A^ = cxd, the above scenario is apparently stable at large A^. 



IV. MASS SPECTRUM AT A^ = 6 

We studied the mass spectrum at A^ = 6, where Eq. (0) predicts the existence of two 
independent mass ratios. In order to extract the other two independent mass values besides 
the fundamental one, we considered the following operators: 

= UabU^d - UadUcb , (26) 

(3) 

Oabcdef = UabUcdUef — UabUcfUgd — UadUcbUef + UadUcfUeb + UafUcbUed — UafUcdUeb , 

(27) 

having respectively the same transformation properties of the two and three particle bound 
states. The mass values M2 and M3 were determined from the large distance behavior of 
the zero space momentum correlation functions constructed with the above operators. In 
practice we found distances d > 1.5 to be large enough to fit the data to the expected 
exponential behavior. In Table and in Fig. [1^ we present the data for the ratios M/Mq, 
M2/M and M^/M, analyzed using the jackknife method. They show good scaling. Fitting 
them to a constant we found 

M/Mg = 0.993(2) , 
M2/M = 1.74(1) , 

M3/M = 2.01(2) . (28) 

This result confirms, within statistical errors of about 1%, the conjectured exact result (^, 
which predicts 

M2/M = 1.73205... , 

Ms/M = 2 . (29) 
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FIGURES 



FIG. 1. Energy versus (3 for = 6. The dashed and dotted hnes represent respectively the 
12^^ order strong coupHng and the third order weak couphng series. 

FIG. 2. Energy versus P for A'^ = 9. The dashed and dotted hnes represent respectively the 
12*^ order strong couphng and the third order weak couphng series. 

FIG. 3. Specific heat and asymptotic scahng test Mq/ Kl^21 for N=3. The dashed line shows 
the analytical prediction (^. Data of C are multiplied by 500. 

FIG. 4. Specific heat and asymptotic scaling test Mq/ Kl^21 for N=6. The dotted line represents 
the IS*'' order strong coupling series of C. The dashed line shows the analytical prediction 
Data of C and the corresponding strong coupling series are multiplied by 500. 

FIG. 5. Specific heat and asymptotic scaling test MqI ^l,21 for N=9. The dotted line represents 
the 13*'' order strong coupling series of C . The dashed line shows the analytical prediction (^). 

FIG. 6. Specific heat and asymptotic scaling test Mg/^l,21 for N=15. The dotted line repre- 
sents the 13*'' order strong coupling series of C. The dashed line shows the analytical prediction 

FIG. 7. Specific heat versus ^g- 
FIG. 8. KL{N,P)/k^2i{N,P) versus /?. 

FIG. 9. Ke{N,(5e)/Ke,2i{N,I3e) versus fdE- 

FIG. 10. M/Mg, M2/M and M3/M versus for iV = 6. Dotted hnes show the exact 
predictions (|) for the ratios M2/M and Ms/M. 
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TABLES 



TABLE L Summary of the immerical results for N=3. 
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1.000(10) 


0.30 


90 


0.43111(2) 


0.176(2) 
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10.41(10) 


0.993(7) 


0.989(9) 


0.315 


120 


0.40400(2) 


0.166(3) 


283(3) 


15.5(3) 


0.981(9) 


0.997(12) 



TABLE IL Summary of the numerical results for N=6. 
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15 
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1.5717(12) 


0.396(12) 
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18 


0.84170(2) 
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0.175 
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0.81060(2) 


0.0399(2) 


2.4806(14) 


0.671(9) 






0.20 


24 


0.77592(2) 


0.0596(3) 


3.0425(18) 


0.816(12) 






0.225 


18 


0.73506(4) 


0.0903(6) 


3.936(3) 








0.25 


18 


0.68234(8) 


0.153(2) 


5.731(9) 


1.332(10) 


0.987(11) 


0.998(15) 


0.26 


24 


0.65537(7) 


0.200(3) 


7.120(11) 


1.560(14) 






0.27 


24 


0.62377(5) 


0.248(2) 


9.407(9) 


1.882(6) 


0.991(7) 


0.998(7) 


0.28 


30 


0.58690(8) 


0.302(4) 


13.68(3) 


2.410(17) 


0.994(6) 


0.997(6) 


0.285 


30 
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0.332(6) 


17.28(3) 


2.836(13) 
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1.007(7) 


0.29 


30 


0.54730(7) 


0.323(4) 


22.31(5) 


3.354(18) 


0.991(8) 


0.992(5) 


0.29 


36 


0.54732(7) 


0.323(5) 


22.38(6) 


3.37(2) 


0.986(6) 


0.996(8) 


0.30 


36 


0.51134(5) 


0.300(4) 


38.27(9) 


4.75(2) 


0.989(4) 


0.996(4) 


0.30 


42 


0.51134(5) 


0.296(3) 


38.30(12) 


4.76(2) 


0.994(6) 
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0.30 


48 


0.51139(8) 


0.288(8) 


38.23(16) 


4.78(5) 


0.999(6) 


0.995(5) 


0.31 


54 


0.48188(4) 


0.256(3) 


65.1(3) 


6.55(5) 


0.986(6) 


1.005(5) 


0.31 


60 


0.48186(4) 


0.258(4) 


65.1(2) 


6.61(4) 


0.992(5) 


1.000(7) 


0.32 


75 


0.45764(3) 


0.227(3) 


108.4(5) 


9.14(7) 


0.996(5) 


0.995(7) 


0.32 


81 


0.45769(3) 


0.230(5) 


107.4(6) 


9.00(11) 


0.997(7) 


1.001(8) 


0.40 


60 


0.34033(5) 












0.50 


60 


0.26335(4) 
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TABLE III. Summary of the numerical results for N=9. 
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Xm 








0.175 


18 


0.812908(11) 


0.03769(13) 


2.4494(12) 


0.672(7) 






0.20 


24 


0.780870(12) 


0.0538(3) 


2.948(12) 


0.792(10) 






0.25 


18 


0.70339(8) 


0.120(2) 


4.868(8) 


1.177(11) 






0.27 


24 


0.65937(13) 


0.205(6) 


6.792(15) 


1.497(17) 


0.980(8) 


0.995(10) 


0.28 


24 


0.62920(12) 


0.268(8) 


8.78(2) 


1.777(13) 


0.990(9) 


0.997(10) 


0.29 


30 


0.58801(13) 


0.410(15) 


13.27(3) 


2.371(14) 


0.994(7) 


0.987(7) 


0.295 


36 


0.56281(12) 


0.457(16) 


17.96(5) 


2.89(2) 


0.994(7) 


0.998(5) 


0.30 


30 


0.53844(21) 


0.420(25) 


25.36(14) 


3.67(3) 






0.30 


36 


0.53845(16) 


0.411(21) 


25.22(12) 


3.63(4) 


0.992(8) 


0.998(8) 


0.30 


42 


0.53847(12) 


0.411(16) 


25.26(9) 


3.68(4) 


0.999(6) 


1.009(6) 


0.31 


42 


0.50028(9) 


0.307(12) 


47.1(2) 


5.39(5) 






0.31 


54 


0.50035(9) 


0.318(12) 


47.2(2) 


5.46(5) 


0.993(9) 


0.995(10) 


0.32 


60 


0.47234(5) 


0.252(10) 


81.6(4) 


7.65(7) 


0.995(4) 


1.006(5) 


TABLE IV. Summary of the numerical results for N=15. 


/? 


L 


E 


C 


Xm 




^g/^w 


^d/^w 


0.20 


24 


0.781395(16) 


0.0515(5) 


2.9392(14) 


0.775(12) 






0.25 


24 


0.70846(3) 


0.100(3) 


4.692(4) 


1.158(14) 






0.28 


24 


0.64990(9) 


0.204(8) 


7.251(11) 


1.561(10) 


0.997(10) 


0.992(9) 


0.29 


24 


0.62134(13) 


0.294(14) 


9.31(2) 


1.857(9) 


0.998(6) 


0.995(7) 


0.30 


30 


0.56809(18) 


0.67(4) 


16.53(5) 


2.750(12) 


0.995(5) 


1.000(6) 


0.31 


45 


0.51202(10) 


0.35(3) 


38.89(16) 


4.78(3) 


0.983(8) 


0.998(7) 








TABLE V. 


Mass spectrum for N=6. 






(3 






M/Mg 




M2/M 




M3/M 


0.29 






0.991(8) 




1.78(3) 




2.02(3) 


0.30 






0.993(3) 




1.74(2) 




2.04(4) 


0.31 






0.992(4) 




1.72(2) 




1.98(3) 


0.32 






0.996(4) 




1.74(2) 




2.06(5) 
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